Jahn- Teller Spectral Fingerprint in Molecular Photoemission: Ceo 
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The hu hole spectral intensity for Ceo — > C^q molecular photoemission is calculated at finite 
temperature by a parameter-free Lanczos diagonalization of the electron-vibration Hamiltonian, 
including the full 8 Hg, 6 Gg, and 2 Ag mode couplings. The computed spectrum at 800 K is in 
striking agreement with gas-phase data. The energy separation of the first main shoulder from the 
main photoemission peak, 230 meV in Ceo, is shown to measure directly and rather generally the 
strength of the final-state Jahn- Teller coupling. 
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Photoemission (PE) from a closed-shell high-symmetry 
molecule, involving Jahn- Teller (JT) effect in the final 
state, is accompanied by characteristic vibronic struc- 
tures in the measured hole spectrum 1] . Although accu- 
rate spectral calculations have recently appeared in the 
chemical literature, e.g. for benzene ^2j, there still is a 
strong need for a qualitative understanding of a more 
general nature, possibly system independent, and infor- 
mative on the nature and the strength of the JT cou- 
pled problem. A recent case in point is that of gas-phase 
Cgo, where an important side peak was reported 230 meV 
above the main PE peak 0,0 (see Fig. As this ex- 
citation energy is in fact larger than all vibrational fre- 
quencies of Ceo (32-197 meV), its nature could not be 
straightforwardly interpreted, leading to a debate . 

In this Letter we present first of all a detailed under- 
standing of the PE structures that are accurately repro- 
duced for Cgo (Fig. through calculations that 
fully include temperature and that contain no adjustable 
parameters. By taking further this calculation to pieces, 
we identify the important ingredients that determine the 
spectrum. The first is symmetry selection: the initial 
orbital symmetry dictated by the suddenly injected bare 
hole extends to the whole spectrum. The second is a 
property of weakly coupled JT vibronic multiplets that 
requires the first excitation energy of the same symmetry 
as the ground state (GS) to rise above the bare vibration 
energy by an amount strictly proportional to the squared 
JT coupling. Specifically in C^g, symmetry requires the 
230 meV peak to correspond to an hu vibronic excita- 
tion above the hu GS. JT transforms the highest Hg 
vibration {hio < 200 meV) of Cgo into an /i„ vibron near 
230 meV. Interestingly, the 20% magnitude ot this up- 
ward shift fingerprints quite accurately the medium-sized 
overall JT couphng of Cgo • As both properties above are 
apparently general, and since the upward shift tracks the 



coupling intensity well beyond the weak-coupling limit, 
this kind of analysis appears of much more general value. 
In particular, the first vibronic excitation energy as seen 
in PE of other weak to medium JT-coupled molecules 
should and does gauge their respective couplings as well. 

We adopt for Ceo the simplest model 0, H, de- 
scribing the JT coupling of the orbitally degenerate /i„ 
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FIG. 1: Experimental PE spectra by Canton et al. Q and by 
Briihwiler et al. compared with the calculation based on 
the ab initio parameters of Ref. (T = 800 K, A'^ticr = 31, 
A'^max ~ 10^, A'sampio = 245). Intensities are normalized to 
unity. Inset: T = 800 K spectra including the single Hg^^ 
mode, and Hg"^^ and IIg^\ All spectra are shifted to move 
the main peak to zero energy and broadened with a 10 meV 
HWHM Gaussian. The errorbar gives the largest statistical 
error introduced by random sampling of the initial states. 



2 



hole with the molecular vibrations. The linear JT cou- 
pling is nonzero for two nondegenerate Ag, six fourfold- 
degenerate Gg and eight fivefold-degenerate Hg vibration 
modes, or 66 vibrations altogether 0, |^ 0| . The result- 
ing h^{A + G + H) JT Hamiltonian thus consists of the 
fivefold hu highest occupied molecular orbital (HOMO), 
of 66 harmonic oscillators, and of the electron-vibration 
(e-v) coupling term 0, 0] : 

H = Hq + Hvih + Hc-v , (1) 

Ho = ShOMO ^ clna'^ma , (2) 
ma 

Hc-v = E f^^9Ajfl^AjC^1^,QtA,,cl^aCm'a- (4) 
o-mm' 

Here mjU label components within the degenerate mul- 
tiplets j counts modes of symmetry A, C^^, are 

Clebsch-Gordan coefficients of the icosahedral group 
Ifi that couple the hu fermion operators cj„^ to a A vibra- 
tion. QiA/i are the dimensionless normal-mode vibration 
coordinates in units of xo{ujAj) — \J^I{'^Kj "t-c) where 
mc is the mass of the C atom), and Pi Ay, the corre- 
sponding conjugate momenta. The additional multiplic- 
ity r = 1,2, needed for Hg vibrations only, labels the two 
separate kinds of Hg coupling allowed under the same 
symmetry [illll). 

In the present calculation we adopt the fully ah ini- 

(r) 

tio numerical values of the e-v coupling parameters 
obtained by a Density Functional (DF) calculation in 
Ref. 11] Hf. We also use the calculated vibrational fre- 
quencies of Ref. H], that are close to the experimental 
frequencies of the Ag and Hg modes, but also cover the 
experimentally unavailable or unreliable Gg modes. E- 
V couplings are dimensionless, each of them normalized 
to its respective vibration quantum, so as to reflect di- 
rectly the relative coupling strength. Numerical factors 
1 i 

^Ag _ 52 _ 52 _ 1 g^j.g i^cludcd to make 

contact with Ref. |6j. 

In the sudden approximation, the angle-integrated PE 
spectrum I{E) is given by Fermi's golden rule 2] 

^(^) = ^E/--H^W' (5) 

mai 

hm.{E) = Y,\{f\cl,\i)\^5{E-Ef+E,) , (6) 
/ 

where \i) represents the starting vibrational eigenstate of 
neutral Cgo, of energy Ei, and |/) represents a vibronic 
JT eigenstate of CgQ, of energy Ef. cl^„\i) is the ini- 
tial bare hole, where the spin-tr electron has been ejected 
from one orbital m of the fivefold degenerate /i„ HOMO, 
but the molecule is still unrelaxed. The factor 1/10 and 



the sum over m and a signifies averaging over the 5 or- 
bital and 2 spin states of the hole. Assuming thermal 
equilibrium we generate the boson numbers {fAj>} of the 
initial neutral molecule by randomly sampling the prob- 
ability distribution P{i) — exp (— i^i/fceT) , where 
Z = Y.i'^^vi-Ei/kBT) and Ei = X]^^Ai^^Aj>- the 
present medium-coupling regime we compute the final- 
state vibronic energies and matrix elements in by 
numerical Lanczos diagonalization of the Hamiltonian 
on the product basis of the five /i„ states times the 
harmonic oscillator ladders. Since these contain an infi- 
nite number of states, some truncation is necessary. As 
we generally want to address high temperatures of the 
molecular beam (T = 800 K used for Cgo, k-oT > 2huj 
of the strongly coupled Hg^"^ mode) a standard trunca- 
tion would involve a far too large Hilbert space size. We 
generate a smarter basis starting from the initial excita- 
tion cJ„crN); and iteratively adding sets of states ("tiers") 
directly coupled to those of the previous tier by matrix 
elements of i?e-vj in a scheme inspired by Ref. [l3|. This 
procedure is iterated iVtier times so that states relevant 
at up to_(2A^tier)-th order in perturbation theory are in- 
cluded [23j . The application of about 350 Lanczos steps 
jl3j | generates a tridiagonal matrix, which provides a well 
converged spectrum. 

The basis for the calculation of Fig. includes about 
720000 states, in particular all the states generated by 
up to three applications of /^c-v to cj-^^\i): the conver- 
gence of the truncated basis is quite good. We repeat the 
whole procedure of basis generation / spectrum calcula- 
tion for each of the A'sampic initial states, and average 
the obtained spectra. The above procedure is carried 
out including all Hg and Gg modes. The Ag modes sep- 
arate out, as they can be solved analytically as a simple 
displaced oscillators, and are included exactly by convo- 
lution. 

Figure H compares the theoretical PE spectrum cal- 
culated at r = 800 K with the measured spectra by 
two separate groups. All experimental features, includ- 
ing the characteristically broad and asymmetric main PE 
peak, the strongest satellite peak around 230 meV, a sec- 
ond weaker peak near 400 meV, and a slow decay up to 
600 meV are remarkably reproduced, based purely on ab 
initio parameters. The model and its ingredients thus 
appear to describe quite accurately the PE spectrum of 
CgQ. We can now take it apart, so as to understand the 
underlying physics to satisfaction. 

We address the two main features of the spectrum, 
namely a) the large broadening and b) the 230 meV 
peak. By eliminating Gg and Ag modes from the cal- 
culation there is no large change in the spectrum. Thus 
only the Hg modes are important. Among Hg modes, 

the largest couplings belong to the low-frequency Hg^'^ 
32 meV quadrupolar radial mode, and to the two high- 
frequency C-C stretch modes Hg'^'^ and Hg^'^ at fku = 181 



3 



and 197 meV). When we include Hg^^ alone we find (in- 
set) that the thermally broadened asymmetric main peak 
is well reproduced. Therefore Hg^^ is the main responsi- 
ble for the large broadening. 

When only Hg^^ is included, all satellite structures 
above 200 meV are absent, and do not reappear even 
if we arbitrarily increase the Hg^"^ coupling to larger and 
larger values. This rules out the possibility to attribute 
the 230 meV structure to an "electronic" splitting 
As it turns out, electronic JT splittings become visible 
in experimental PE spectra, such as those of Fe(C0)5 
[l^ . and in theoretical spectra, like those of Martinelli 
et al. |l5l | only for very strong e-v coupling, whereas in 
Ceo coupling is intermediate at most. That suggests that 
the 230 meV peak in Ceo should rather be associated to 
the high-frequency modes Hg'^^ and Hg^\ However, the 
connection between their frequencies < 200 meV and the 
230 meV peak position is far from obvious. To clarify 

(7) 

this, we calculate the spectrum by including Hg ' and 
alone (inset of Fig. [TJ. The result is a much nar- 
rower spectrum (i?g^^ is now omitted), showing a high- 
energy satellite at a significantly (10%) higher energy 
than huJHgS- This increase in vibronic excitation energy 
above the purely vibrational hio for the u=l-phonon ex- 
citation of a trivial displaced oscillator is in fact a char- 
acteristic signature of the dynamic JT effect, and can be 

(7) (S) 

understood on the basis of symmetry. For Hg and Hg , 
fko is almost three times the thermal energy: the main 
initial-state contribution to the spectrum comes from |i) 
= the neutral molecule GS, of symmetry Ag. The bare 
hole cJ^ctN) therefore carries the orbital /i„ symmetry of 
the HOMO. This state has nonzero overlap strictly only 
with final states |/) of the same symmetry hu- This single 
symmetry channel thus represents the only contributor to 
the low-temperature PE spectrum according to JHl 24] . 
The blue shift of the vibronic peak relative to the bare 
vibration is now due to two general properties (in fact 
common to all dynamical JT systems): (i) the coupling 
independence of the average excitation energy in any vi- 
bronic multiplet at weak coupling (the leading change is 
of order g"^) [ig; (ii) the increasing "repulsion" that the 
GS level exerts for increasing g on all excited states of 
the same symmetry, which pushes them upward. 

Both concepts are illustrated in Fig. [3 depicting the vi- 
bronic excitation energies (with special attention to the 
hu states) for the simplest case of a single Hg mode, 
as a function of the e-v coupling strength g^. The ap- 
proximate coupling independence of the average excita- 
tion energy of the multiplet of states derived from the 
i;=l-phonon manifold is apparent. Within this multiplet 
(made of huXHg = au + tiu+t2u + '2'gu + '2hu), the /i„ vi- 
bronic states are the sole that have another state of the 
same symmetry (the GS itself) lower in energy "push- 
ing" them upward. As a result, for small but increasing 
coupling, the hu excited vibronic states, move necessarily 
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FIG. 2: Vibronic excitation energies for a single Hg mode 
coupled to the hu level, as a function of coupling, for a ra- 
tio g'^Vff''^'' characteristic of the H'^^^^ modes. The average 
refers to the multiplet originating from the 1-phonon states. 
Bold: hu states, the only visible states in T = PE. 



upward in energy above hwH^ ■ Moreover the amount of 
shift is at weak coupling proportional to g^ and reflects 
the strength of the coupling. None of these results seems 
specific to the h^H JT model of Ceo. For example, iden- 
tical conclusions apply to the e®E problem J6]. Indeed 
a 15% blue shift of the first vibronic satellite is clearly 
observable in the PE spectra of benzene fl7^. 

In Ceo we must consider the two high-frequency modes 
H^P and ijf a.t huj ^ 181 and 197 meV). The couplings 
for these two modes alone accounts for a 10%, or roughly 

(7 8~1 

20 meV shift of the main H^ ' vibronic satellite above 
hjjj. This is about half of the total observed and calcu- 
lated shift of ~ 40 meV of the 230 meV satellite relative 
to hw. The remainder is a collective effect of all other 
modes simultaneously interacting with the HOMO, and 
pushing their respective v = \ satellites (weak and thus 
invisible at large temperature) upward. All goes as if ef- 
fectively the coupling of the h'^^^^ modes were larger, 
i.e. effectively moving to the right of the vertical dashed 
line in Fig. |5| 

We have thus attained the following understanding 
of the 800 K PE spectrum of Ceo: (i) the strong 
230 meV satellite (and the weak 400 meV too) is due to 

(7 /Si 

/i„-symmetry vibrons derived from Hg modes, blue 
shifted by repulsive coupling to the /i„ GS; (ii) the large 
spectral broadening is due mainly to h]^\ Clearly, this 
broadening of the 800 K spectrum prevents direct ex- 
traction of further detailed information about the actual 
e-v couplings of the fuUerene cation. Our calculations, 
when repeated for lower temperatures, indicate a wealth 
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FIG. 3: Temperature dependence of the PE, computed in- 
cluding all modes of Ceo for T = 0, T = 300 K and T = 800 K. 
Labels A'^',, indicate the leading component in the vibronic 
states. Inset: computed spectra (T = 800 K) based on the 
same DF coupling parameters divided by \/2, original, 
and multiplied by \/2, compared to experiment 0). 



of structures that should emerge. Figure 13 shows that al- 
ready a measurement carried out at 300 K would provide 
enough detail to get precise indications about the cou- 
plings of individual modes. For example, the intensity 
ratios between the main peak and the subsequent peaks 
at 35 and 65 meV is a direct measure of gHgi- A cool 
molecular beam experiment would be extremely use- 
ful. Also in view of future cool-beam data it is meaning- 
ful to address the question of how accurate the calculated 
DF e-v couplings used here really are. For tiu electrons 
coupled to Hg vibrations, the calculated DF couplings 
were argued for example to be too weak by a substantial 
amount, by comparison to the larger couplings needed in 
order to fit Cqq PE In our HOMO case, the excel- 
lent agreement of calculations with experiment suggests 
that there is no such failure. Figure OK inset) shows sev- 
eral PE calculations with rigidly rescaled e-v parameters. 
Comparison with experiment indicates that the DF cou- 
plings are, on the whole, quite close (if slightly on the 
weak side) to the actual couplings of C^q. Much larger 
JT coupling values like those obtained using different ap- 
proaches ,20| are ruled out. That is very good news, for 
it gives further confidence in the future possibility to rely 
solidly on DF for a good description of the intra-fuUerene 
chemical bond, and of its perturbations. 

In conclusion, we present here a detailed calculation 
of the PE spectrum of a high-symmetry molecule, Cgo, 



which is JT active in its one-hole state. We find excellent 
agreement with experiment, confirming for the first time 
the accuracy of the e-v couplings as extracted from DF 
calculations. Perhaps more importantly, we reach a fresh 
understanding of the PE spectrum, with its width, its 
vibronic peak, and the reason for its blue shift relative 
to the bare vibrations. The blue shift of the one-vibron 
satellite is argued to represent a more general characteris- 
tic fingerprint of JT systems in the weak-to-intermediate 
coupling regime, and is found for example in the e (E) E 
JT of benzene as well. 

We are indebted to G.P. Brivio, P. Briihwiler, A. Del 
Monte, D. Galli, M. Lucders, L. Molinari, G. Onida, 
F. Parmigiani, A. Parola, and G.E. Santoro, for use- 
ful discussion. This work was supported by the Eu- 
ropean Union, contracts ERBFMRXCT970155 (TMR 
FULPROP) as well as by MIUR COFINOl, and FIRB 
RBAU017S8R operated by INFM. 
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